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FOREWORD 


The  Naval  Surface  Weapons  Center's  Independent  Research  (IR)  Program 
is  designed  to  stimulate  original  work  of  a  basic  nature  and  increase 
competence  in  all  fields  of  science  and  technology  relevant  to  the 
Center's  mission. 

This  report  documents  one  facet  of  the  engineering  studies  performed 
by  the  author  under  an  IR  project  titled,  "Adaptive  Signal  Detection," 
during  the  first  half  of  1983. 
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SUMMARY 


The  problem  of  weak  signal  reception  in  non-Gaussian  noise  has  been 
shown  to  have  an  optimal  solution  the  structure  of  which  involves  a 
no-memory  nonlinearity  followed  by  a  matched  filter.  The  shape  of  the 
nonlinearity  is  completely  determined  by  the  noise  amplitude 
probability  density  f unction^(pdf ) The  assumptions  underlying  this 
theory  are  examined  and  applications  to  radio  communication  in 
atmospheric  noise  are  considered.  Since  the  likelihood  ratio  test  on 
which  the  theory  is  based  is  not  a  robust  procedure,  the  receiver 
system  must  adjust  the  shape  of  the  nonlinearity  to  account  for  the 
statistics  of  the  prevailing  noise  environment.  ^ 

INTRODUCTION 

The  fidelity  of  communications  between  a  transmitter  and  receiver  at 
opposite  ends  of  an  additive  noise  channel  can  be  improved  (1)  by 
increasing  the  ratio  of  signal  power  to  noise  power,  (2)  by  changing 
the  form  of  the  signal  while  holding  power  constant,  or  (3)  by 
designing  better  noise  immunity  into  the  receiver.  This  publication 
concerns  the  third  option  with  emphasis  on  the  case  in  which  the 
received  signal-to-noise  ratio  (SNR)  is  low. 

Kotel'nikov  established  four  categories  for  the  classification  of 
noises  that  may  interfere  with  radio  communications:  (A)  sinusoidal 
noise,  (B)  impulse  noise,  (C)  normal  noise,  and  (D)  super-positions  of 
B-  and  C-type  processes.^  Impulse  noise  is  characterized  by 
intermittent  pulses  of  mean  duration  much  less  than  the  mean 
interarrival  time.  Normal  noise  can  be  understood  as  a  sum  of 
independent  impulsive  processes  so  numerous  that  the  mean  interarrival 
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time  for  pulses  in  the  sum  process  is  much  less  than  pulse  duration. 
Then  the  instantaneous  amplitude  of  the  sum  process  tends  to  be 
normally  distributed  (so  long  as  the  variance  of  the  pulse  height  is 
not  too  great) .  Type-D  noises  are  those  in  which  the  number  of 
independent  impulsive  processes  is  not  sufficient,  or  the  mean  rate  of 
pulsing  is  too  slow,  for  the  central  limit  theorems  to  apply.  Thus 
the  successive  sample  values  of  the  type-D  process  may  be  independent 
and  identically  distributed  according  to  a  law  different  from  the 
normal.  Although  the  four  categories  of  Hotel' nikov  do  not  exhaust 
the  range  of  possibilities,  they  serve  to  motivate  a  fairly  general 
treatment  of  the  problem  of  designing  a  communications  receiver  with 
optimal  noise  immunity.  It  should  be  noted  that  a  second  dimension  of 
classification  is  always  implicit  in  the  characterization  of 
interference  sources,  as  they  are  either  environmental  or  synthetic 
(system-generated) . 

An  understanding  of  environmental  noise  in  radio  communications  at 

frequencies  up  through  several  Megahertz  evolved  slowly.  Eccles 

(1912)  is  cited  by  Austin  as  being  tirst  to  propose  that  broadband 

atmospheric  interference  originated  principally  in  lightning  strokes, 

2 

both  locally  and  at  great  distances.  DeGroot  (1917)  suggested  that 
the  magnetosphere,  or  exosphere,  is  the  medium  through  which  solar 
winds  generate  atmospherics.  Both  hypotheses  are  accepted  in  current 
theory.  The  first  demonstration  that  "whistlers"  propagate  in  the 
exosphere  was  made  in  1959  using  Vanguard  III.^  Shepelavey  modelled 
atmospheric  noise  up  to  20  MHz  as  the  sum  of  impulsive  and  normal 

4 

processes.  In  other  words,  if  X  is  the  amplitude  due  to 
atmospheric  noise,  one  has  X  *  Y  +  aZ,  where  Z  is  standard  normal,  o 
is  the  r.m.s.  normal  fluctuation,  and  Y  is  the  impulsive  contribution 
described  by  the  two-sided  power-Ray leigh  distribution. 


P ( I Yl >y)  *  exp-(by)a,  y>0 


(1)  , 


p  *  fraction  of  time  x  is  exceeded 


FIGURE  1.  EXCEEDENCE  CURVE 


facilitated  by  plotting  the  coordinates 


u  *  -log(-lnp) 


(on  the  horizontal  scale)  and 
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(on  the  vertical  scale),  as  shown  in  Figure  1,  where 


p  *  P(  lXi>x) 


is  the  probability  that  the  absolute  sample  value  exceeds  x.  The 
result  is  known  as  an  exceedence  curve  and  is  illustrated  above. 

Field  and  Lewinstein  showed  that,  after  setting  asl/4  and  adjusting 
b  to  give  the  right  division  of  mean  power  among  the  two  component 
processes,  the  convolution  of  the  power-Rayleigh  and  normal  densities 
gives  a  density,  and  hence  an  exceedence  curve,  which  matches  the  ELF 
and  VLF  data  extremely  well.^ 

Impulse  noise  may  also  be  synthetic.  If  communications  receivers  are 
placed  near  radars,  impulsive  interferences  may  occur  due  to  direct 
pickup,  feedback  through  common  power  supplies,  or  shock  excitation  of 
tuned  circuits  of  the  receiver  by  large  modulator-generated  d.c. 
pulses.**  "Ignition  noise"  refers  to  one  of  the  earlier-known  types 
of  synthetic  impulsive  interference. 

The  treatment  prescribed  for  impulse  noise  by  a  number  of  authors  in 
the  1930's  and  -40' s  was  to  "silence"  or  "blank"  the  receiver.  These 
and  other  terms  symbolized  steps  in  the  development  of  the  limiter 
concept.  A  short  duration  pulse  at  the  input  of  a  tuned  circuit  in  a 
narrowband  receiver  system  will  produce  a  damped  sinusoidal  response 
with  damping  rate  essentially  determined  by  the  "Q”  of  the  circuit. 

The  communication  of  data  will  be  interrupted  until  the  "ring  down" 
waveform  subsides.  If  a  broadband  amplifier  feeding  the  tuned  circuit 
is  made  to  saturate  at  a  level  which  is  intermediate  between  the 
height  of  the  input  pulse  and  the  nominal  signal  amplitude,  then  the 
effect  of  the  noise  pulse  is  mitigated  while  the  signal  passes 
undistorted.  Other  limiters  are  possible.  A  "hole-puncher"  is 
obtained  bi  makin'  .he  amplifier  output  drop  to  zero  when  the  limiting 
level  is  exc.-J«d.  A  limiter  would  also  be  useful  in  a  broadband 
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receiver  system  which  integrates  the  input  over  a  fixed  interval  to 
discriminate  between  positive  and  negative  signal  pulses  of  level  +S. 
Brief  noise  pulses  of  amplitude  +A,  A>>S,  will  strongly  influence 
the  result  of  the  integration.  The  suppression  of  large  excursions 
before  integration  improves  fidelity.  In  general,  a  limiter  is  any 
memoryless  device  which,  given  an  input  x,  immediately  outputs  g(x),  a 
function  that  is  strictly  bounded  above  and  below.  If  a  limiter  cuts 
off  at  values  of  x  which  lie  in  the  range  of  probable  signal 
amplitudes,  then  the  spectrum  of  g(x)  is  related  to  that  of  x  in  a 
generally  complicated  way.  The  term  detector  usually  refers  to  a 
no-memory  nonlinearity  which  is  inserted  in  the  receiver  to  effect 
desired  changes  in  the  spectrum.  For  example,  a  diode  or  square  law 
detector  is  used  to  demodulate  an  AM  signal.  It  would  not  be  wholly 
inconsistent  with  current  usage  to  say  that  a  generalized  limiter  is  a 
detector  which  suppresses  noise  at  some  acceptable  cost  in  signal 
distortion.  Middleton  gives  detailed  accounts  of  how  different  types 
of  detectors  perform  in  AM  and  FM  receivers  across  the  range  of 
SNR ' s . 7 


LIKELIHOOD  RATIO  DETECTORS 


Until  the  mid-1960's,  the  treatment  of  detectors  for  noise  immunity 
enhancement  was  ad  hoc.  Development  of  a  general  theory  based  on  the 
optimal  properties  of  likelihood  ratio  tests  appears  to  begin  with 
Rappaport  and  Kurz  who  showed  that,  for  binary  synchronous 
communication  in  non-Gaussian  noise,  the  optimal  receiver  consists  of 
a  two-input  nonlinear  device  followed  by  an  integrator  and  a  decision 

O 

box.  Consider  the  problem  of  binary  synchronous  communication  in 
the  presence  of  additive  noise  y(t)  having  pdf  f(y).  It  is  presumed 
that  the  transmitted  signal  is  s^(t),  i  =  1  or  2,  -T/2  <  t  <  T/2.  The 
objective  is  to  design  a  receiver  which,  based  on  the  sampled  data  {xn 
■  x(tn),  n*0,l,...,  N> ,  where 


xn  *  sin  +  yn 


(2)  , 


5 


Vo'. 


V>i\  <\ 
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will  decide  correctly  between  the  two  hypotheses.  The  likelihood 
ratio  can  be  expressed  as  the  product 


A 


N 


-  L 


t(xn  ~  sln> 
f(xn  "  s2n) 


(3) 


of  component  ratios  under  certain  conditions  which  are:  either  (i) 
the  type-D  or  -B  random  process  described  by  f(x)  has  a  mean 
interarrival  time  much  greater  than  T,  or  (ii)  the  system  bandwidth 
greatly  exceeds  the  signal  bandwidth.  These  two  conditions  serve  to 
justify  the  factorization  of  the  joint  density  which,  in  the  strict 
sense,  follows  from  the  independence  of  the  y  .  Defining 


Jn 


f  (x  - 
-  log  n 


sln} 


f(xn  -  s2n) 


(4) 


one  has 


logA  = 


E 

n=0 


Jn 


(5) 


7 

which  reduces  the  likelihood  ratio  test  A  >  C  to 


N  ? 

J  Ln> 
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The  conditions  which  justified  factorization  of  the  joint  density  of 
the  N+l  sample  values  now  allow  a  normal  approximation  to  the  distri 
bution  of  logA ,  since  every  Ln  is  well-behaved  for  continuous 
densities  f.  The  proposed  receiver  is  block  diagrammed  in  Figure  2. 

For  square  antipodal  signals. 


s,  (t) 


Sn(  t  ) 

2 x  * 

FIGURE  2.  LIKELIHOOD  RATIO  RECEIVER 
( S ,  -T/2  <  t  <  T/2 


si  (t) 


!S,  -T/2  <  t  < 
0  otherwise; 


s2  (t)  =  -sx  (t)  (7)  , 

the  threshold  C-l  corresponds  to  selection  of  s^  when  the  likelihood 
ratio  statistic  is  positive,  s2  when  negative.  The  error  probability 
is 


Pe  =  P  (logA<,0  '  i=l)  =  (1/2) erfc  (p/VT) 
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where 


p  =  E  (logA  |  i=l)  /  (Var  (logA ) ) 1//2  (9) 

is  the  ratio  of  the  mean  to  the  standard  deviation  at  the  integrator 
output,  which  can  be  called  the  output  SNR.  Note  that  logA  has  the 
same  variance  for  i  ■  1  or  2  when  (7)  applies. 

Antonov  (1966)  proposed  a  more  interesting  treatment  of  the  problem, 
beginning  again  with  equation  (3) ,  but  now  testing  the  hypothesis 
s(t)=0  versus  the  simple  alternative  in  which  s(t)  is  any  nonzero  real 

o 

waveform.  The  first  sum  in  the  equation 


N  N 

logA  =  T  logf(xn-sn)  -  Y*  logf(xn) 

n=0  £o  dO) 


is  expanded  in  a  Taylor  series  about  the  points  sn*  Interchanging 
the  order  of  summation  and  collecting  terms,  one  has 


logA  =  V  il^-ll  sK  -ill  logf(y) 

k-1  kJ  n-0  n  iyk 

A  00 

*  2  <Pk 
k=l 


(ID  . 


For  a  weak  signal, 


logA  sng(xn) 


I  ■»  •**  V 
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to  first  order,  where 


g  (x)  ■  -  SL  logf  (x) 
dx 


(13) 


defines  a  no-memory  nonlinearity  (a  "lagless  fourpole") .  The 
operation  of  the  optimal  first  order  receiver  therefore  consists  in 
(1)  passing  xn  through  a  detector  to  produce  9(xn)  and  (2)  summing 
the  point-by-point  products  of  g(xn)  and  the  stored  values  of  the 
known  signal.  The  second  operation  describes  the  response  of  a 
matched  filter.  Not  accidentally,  if  f (x)  is  the  standard  normal  pdf, 
equation  (13)  reduces  to  9 (xn)  *  xn»  and  the  optimal  first  order 
receiver  computes  the  sum  of  terms  xnsn,  i.e.,  the  solution  for 
white  Gaussian  noise  is  just  a  matched  filter. 


FIGURE  3.  ANTONOV'S  OPTIMUM  DETECTION  RECEIVER 

Figure  3  shows  Antonov's  optimum  detection  receiver  for  a  completely 
known  weak  signal,  where  INF  is  the  nonlinear  device,  OLF  is  the 
optimum  linear  filter,  and  DD  compares  9^  »  <p  to  a  threshold  9*  to 
decide  between  HQ  (noise  only)  and  (signal  present) .  It  remains  to 
show  that  the  operating  characteristic  of  this  receiver  is  in  fact 
better  than  that  of  a  matched  filter  (OLF)  alone.  To  this  end  let  ^ 
denote  the  matched  filter  output  under  the  tacit  assumption  that  INF 
is  deleted  [g(x)=x],  Then  will  be  compared  to  a  threshold  ty*.  When 
N  is  sufficiently  large,  both  and  9  are  approximately  normal  under 
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both  Hq  and  H^;  an  solution  requires  computation  of  the 
conditional  means  iriances.  Let  EQ  =  E(»  IHQ),  E^  *  E(«  ih^) 
similarly  for  the  ice  operator.  Three  assumptions  make  the 
computation  relati;traightforward.  (1)  The  density  f(y)  is 
even  function;  i.eJ'ff-y).  (2)  The  signal  autocorrelation 
is  much  less  than  iterval  between  samples  so  that 


Z2  A 

sn)  sn  =  S . 


(3)  The  signal  is  .nough  that 


g(sn+xnxn)  +  sng '  (xn) 


is  generally  valid 


These  assumptions  :o  the  results 


( E°k<P  ( 

\Elkcp  \V 


kSE0g*  s 


and 


/varDkcp  A  .  /k2SEQg2  a2 


\Varik(p  4 ,)  \k2SE0g2  c2 


i  and 

an 

width 


(14) 


(15)  , 


10 
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where  k  is  a  constant  that  gives  k«p  the  same  units  as  i|> .  Of 
particular  interest  is  the  result  for  the  mean  of  9  under  H^, 
derived  with  the  aid  of  (3) .  The  key  quantity  is 


E092 


00 

J  f(x)g2(x)dx  *  E0<p2/s 
—00 


(16) 


Defining 


R  =  S/a2 


(17) 


the  SNR  at  the  output  of  the  matched  filter  when  ^  is  considered, 
simple  algebra  shows  that 


Eo<P 


2  = 


Rn 


where  the  number 


n  *  o2 


00 

/  £(x,[l;l09f(x)] 


dx 


(18) 


(19) 


measures  the  effective  SNR  improvement  gained  by  using  INF.^®  Using 
the  Schwartz  inequality,  it  can  be  shown  that 


NSWC  TR  83-236 


with  equality  if  and  only  if  the  noise  is  Gaussian.  The  false  alarm 
and  detection  probabilities  are  given  explicitly  by 


Q0  "  p  (9  *  f*  IHo^  =  1  ”  *(<p*/VnS) 


(21 


and 


Qd  =  p  (<P  >  9**^1)  =  1  ■  <t  (cp */v^R  -  VnR) 


(22 


respectively,  where 


*  ( z ) 


(2w)’1/2exp(-t2/2)dt. 


Since  <p*  and  are  both  adjusted  to  give  the  same  Qq  in  the 
Neyman-Pearson  test,  one  has 


(23 


Now  the  detection  probabilities  are  given  by  (22)  and 


Qd  *  P(\|)><j>*/V?r  Hl) 

«  1  -  *  (<p  Vn?R  -  VS) 


(24 
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The  assertion  reduces  to 

q>*  <p*/n  -  (l-l/VrT)R  (25) 

completing  the  argument. 


ADAPTATION  TO  CHANGING  ENVIRONMENTS 

The  optimality  theory  outlined  in  the  preceding  paragraphs  presumed 
foreknowledge  of  the  first  order  density  f (y) .  This  knowledge  would 
allow  the  designer  to  construct  a  detector  with  the  appropriate  g (x) 
-dlogf (x)/dx.  A  question  arises  concerning  the  effect  of  temporal 
variation  in  the  density  f(y;t),  where  t  is  any  time  in  the  life  cycl 
of  the  receiver  system.  This  question  concerns  the  robustness  of  the 
likelihood  ratio  test  on  which  the  theory  is  based/  as  a  statistical 
test  or  estimation  procedure  is  said  to  be  robust  if  its  accuracy  is 
not  too  sensitive  to  flaws  in  the  underlying  stochastic  model.  No 
satisfactory  answer  of  broad  generality  is  apparent.  Consequently, 
the  designer  must  address  the  changing  environment  of  the  receiver  by 
designing  subsystems  which  make  the  receiver  adapt  to  the  noise.  In 
other  words,  the  optimum  nonlinear  receiver  is  inherently  an  adaptive 
device  which  adjusts  g(x)  to  reflect  the  recent  history  of  the  noise. 

In  order  for  this  solution  to  be  workable,  a  degree  of  short-term 
stability  is  required  of  the  noise.  To  illustrate  this  point, 
consider  the  case  of  a  multimode  ergodic  process  {Yfc,  -«><t<<*>} 
for  which 
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for  every  t.  The  index  j  identifies  the  state  of  the  process  and,  for 
J  sufficiently  large,  the  postulate  (26)  will  almost  surely  be  valid 
modelling  assumption.  If  f(y;t)  belongs  to  a  family  of  densities 
£(y;Q)  parameterized  by  a  vector  j),  then  (26)  amounts  to  saying 
that  the  (continuous)  space  containing  every  possible  j)  can  be 
partitioned  into  J  disjoint  regions  each  small  enough  that  variation 
within  a  single  region  is  of  negligible  impact  on  the  density.  Then 
Q(t)  is  approximated  by  the  sequence  of  discrete  states  ljt>, 
and  f(y;t)  =  f(y;JJ(t))  is  approximated  by  f  j  (y)  .  If  the  mean  time 
between  transitions  from  one  state  to  the  next  generally  exceeds  the 
maximum  length  of  a  message,  the  detector  which  uses  g^  (x)  =  -dlogf(x)/ 
dx  will  be  suboptimal  only  on  rare  occasions. 

Given  a  parametric  pdf  which  evolves  slowly  enough  to  permit 
application  of  the  theory,  the  designer  needs  to  provide  receiver 
subsystems  which  can  estimate  the  distribution  parameter  and  use  the 
estimate  to  modify  the  detector  as  required.  A  broad  range  of 
procedures  apply  to  the  problem  of  probability  density  estimation.  It 
is  assumed  that  the  channel  is  unused  most  of  the  time,  so  that  by 
estimating  the  present  density  from  the  recent  sample  values  of  the 
input,  the  receiver  really  does  obtain  f  (y) ,  and  not  the  density  of 
y(t)+s(t).  Since  s(t)  is  known,  a  corresponding  pdf. 


fs(x)  =  P (s ( t) <x) 


(27)  , 


is  known;  and  the  density  f_.„(x)  of  the  signal  plus  noise  can  in 

S*ry 

principle  be  used  to  get  f(y)  by  solving  the  integral  equation 


00 

fs+y(x)  *  f  fs  (x-y)  f  (y)dy 
—00 


(28) 
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for  f(y).  Yet  if  the  receiver  does  not  have  access  to  a  fairly 
sophisticated  computer/  one  would  not  want  to  address  this  issue. 
Besides/  if  message  length  is  much  less  than  the  period  of  time  over 
which  the  input  is  sampled  to  estimate  j),  the  presence  or  absence  of 
the  weak  signal  in  some  of  the  samples  will  not  strongly  influence  the 
result. 

Three  approaches  to  pdf  estimation  were  considered  by  Evans,  who 
addressed  the  present  problem  of  finding  g(x)  for  Antonov's 
receiver. ^  Distribution— Independent  Linear  Smoothing  computes 


+u/2 

f  (x)  =  f  wn(x-y)dFn(y) 

J-  u/2 

from  the  accumulated  histogram  Fn  of  n  sample  values.  The  smoothing 
function  wn  is  suitably  well  behaved  and  even  on  the  interval  (-u/2, 
+u/2)  and  it  approaches  a  delta  function  as  n+®.  Partially 
Distribution-Independent  Linear  Smoothing  uses  a  function  wn  in  (29) 
which  has  been  selected  to  minimize  the  variance  of  fn  over  a  given 
family  of  density  functions. 

Evans  used  the  term  Distribution-Dependent  Smoothing  to  describe  the 
case  in  which  one  assumes  a  model  f(x;J3)  and  estimates  using 
techniques  such  as  maximum  likelihood  or  the  method  of  moments. 
Application  of  these  two  techniques  to  the  power-Rayleigh  distribution 
introduced  in  Equation  (1)  is  a  straightforward  exercise.*  If  Y  is 
the  random  variable  and  a*l/m,  the  first  two  moments  are  EY  =  ml/b 
and  EY^  *  b“^r (2m+l) .  The  method  of  moments  prescribes  the 
estimates  &  =  ni!/Y,  where  Y  is  the  sample  mean,  and  m  such  that 


*Here  Y  corresponds  to  the  lYl  of  equation  (1) . 
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Y2/  ( Y)  2  =  r  (2m+l) /  (nil )  2 


(30)  , 


■"*2  A 

where  Y  is  the  sample  second  moment.  Solving  (30)  for  m,  the  result 

can  be  substituted  into  the  equation  for  £. 

The  maximum  likelihood  estimators  of  a  and  b  are  considerably  more 
complicated.  One  ends  up  with  the  simultaneous  equations 


*  A 

n A  -  £  (1  -  (ftYi)a)lnYi 

i=l 

and 

n  -  £  [fivj2 
i*l 


for  sample  values  Y^,...#Y  . 

Without  an  assumed  parametric  model  and  beginning  with  a  smoothed/ 
normalized  histogram  p(x)  that  has  been  derived  from  old  data/  an 
updated  pdf  estimate  q(x)  could  be  obtained  using  the  method  of 
minimum  cross-entropy  described  by  Shore.  Let  the  true  pdf  be 
f (x) .  A  set  of  linear  equality  constraints/ 


J f  (x)  e  j  (x) 


ax 


(31)  / 
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together  with 

Jq  (x)dx  *  1  (32) , 

lead  to  the  result 


K 

q(x)  ■  p (x) exp ( -X  -  flje j  (x)  ) 

j=l 


(33)  , 


for  unknown  Lagrange  multipliers  X  and  .  The  method  of 

minimum  cross-entropy  selects  that  q(x)  which  minimizes 


H  (q/P)  =  ^q(x)  log|q(x)/p(x)J  dx 


(34) 


The  cross-entropy  at  the  minimum  can  be  expressed  in  terms  of  the 
Lagrange  multipliers  as 


K 

minH(q,p)  *  -X  -  ^  Bjdj 

j=l 


(35)  . 


For  example/  one  could  have/  for  K=2f  that  d^  and  d2  are  the  first 
two  moments  of  the  distribution.  Then  the  equations 


Jp(x)  (x-di)exp(-B;]X 
/p(x)  (x2-di)exp(-fiix 


fl2x2)dx  =  0/ 

-  fl2x2)dx  =  Of 
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and 


X  *  log  y*p (x) exp (-B^x  -  13 2X 2 ) x 


could  be  solved  for  the  Lagrange  multipliers.  Unfortunately,  this  too 
is  a  rather  computation-intensive  procedure. 

CONCLUDING  REMARKS 

The  optimality  theory  outlined  above  can  be  illustrated  with  a 
reasonably  simple  example.  Referring  back  to  Shepelavey's  atmospheric 
noise  model  (page  2) ,  let  the  sample  value  of  the  noise  be  X  =  Y  +  cZ 
with  the  two-sided  power-Ray leigh  distribution  of  Y  having  a=l« 

(Let  c  denote  the  r.m.s.  normal  noise.)  Then 


f (y)  -  (l/2)bexp(-biyi )  (36) 

2  -2 

is  the  bilateral  exponential  density  with  EY  *  2b  .  For  c=0, 
the  optimal  detector  has 

g ( 0 ) (x)  «  bsgn(x)  (37), 

which  describes  a  "hard  limiter";  and,  for  c*l, 

g (1) (x)  ■  x 

as  noted  previously.  Every  intermediate  case 

g  (x)  ,  0<c<  1, 
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requires  evaluation  of  a  convolution  integral 


fx<x) 


00 

/2  2 

f(x-y)e~y  /2c  dy/V^iTc 

—00 


(38)  . 


It  is  a  straightforward  exercise  to  program  a  desktop  computer  to 
evaluate  (38)  for  a  number  of  cases  subject  to  the  mean  power 
constraint 


c2  +  2/b2  =  1 


(39)  . 


(c )  2 

Figure  4  shows  the  shape  of  the  optimal  g  (x)  for  c  =  0.975, 

0.50,  0.25,  0.125,  and  0.05.  These  graphs  appear  to  be  consistent 

with  the  requirement  that  g(x)  approach  (37)  as  c+0.  The  table  of 
2  2 

SIGMA  ■  c  versus  IMP  *  n,  inset  in  the  figure,  gives  the 

computed  SNR  improvement  factors  obtained  from  the  integral  in 

2 

equation  (19) ,  with  o  *1.  The  accuracy  of  the  computation  is 
clearly  not  too  good  in  this  respect.  (For  c=l  the  result  should  be 
r|*l.)  Yet  the  trend  is  towards  increasing  n  as  the  non-Gaussian 
fraction  increases.  For  b  =  VT,  the  exact  solution  is  n=2. 

Using  smoothed  sample  pdf's  derived  from  ELF  atmospheric  noise  data, 
Evans  and  Griffiths  found  the  SNR  improvement  factor  to  be  generally 
on  the  order  of  a  few  dB.  In  light  of  the  difficulties  involved  in 
designing  an  adaptive  optimal  detector,  they  chose  to  consider  various 
kinds  of  fixed  limiters  for  their  hypothetical  ELF  receiver.  The 
rapid  advances  in  microcomputer  technology  over  the  intervening 
decade,  however,  could  justify  a  re-evaluation  of  their  conservative 
conclusion. 
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